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Cn , We discuss the importance of observing supernova neutrinos. By analyzing the SN1987A 

OO 1 observations of Kamiokande-II, 1MB and Baksan, we show that they provide a 2.5(t 

^D ' support to the standard scenario for the explosion. We discuss in this context the use of 

CO ' neutrinos as trigger for the search of the gravity wave impulsive emission. We derive a 

bound on the neutrino mass using the SN1987A data and argue, using simulated data, 
that a future galactic supernova could probe the sub-eV region. 
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Introduction Core collapse supernovae (type II, lb and Ic) occur when the pro- 
genitor has a mass > 8 Mq and originate compact remnants: neutron stars, black 
holes and possibly hybrid (=quark core) stars. The formation of such an object re- 
quires to carry away a huge binding energy, several times 10^'^ erg. It is well known 
that the role of carriers is played mainly by neutrinos; but more importantl y for 
us, and according to the standard scenario of core collapse supernova explosioinJ^, 
neutrinos play also a fundamental role in driving the explosion. They deposit energy 
that can revive the shock, which will eventually cause the expulsion of the external 
layers of the star. 

The current scenario of neutrino emission is based on two main phases of neu- 
trino emission. The first one, called accretion phase, entails 10-20% of the total 
energy. It is characterized by a very high neutrino luminosity and is directly related 
to the matter which is accreted over the proto-neutron star through the stalled su- 
pernova shock wave. The other phase is called cooling phase; the neutrinos escape 
slowly from the proto-neutron star, releasing the remaining 80-90% of the energy. 
Only two analyses of SN1987A data included both emission phases: the analysis of 
2001 by Loredo and Lamb^and the recent one due to our group*^. The most relevant 
modification of this last analysis is the improvement of the model of D^ emission 
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Fig. 1. Directions of events events in a Cherenkov detector of 32 kton: the ~ 300 elastic scattering 
events are visible over the ~ 5000 inverse beta decay events, expected from the galactic center. [We 
use the Lambert projection: the points of the unit sphere— i.e., the possible directions-identified 
by n = {sgc^,SQS^,cg), are mapped into the circle of radius 2, whose points are identified by 
(u,v) = (c0, s^)y'2(l — eg); namely, u = n^^ljiX — n^) and v = Uy^ljiX — n^,). The Lambert 
projection conserves the areas, dO = d(j)dcg = d?a = dudv.] 

that we are going to describe in some detail in the following^ Remarkably, both 
these analyses claim an evidence of the phase of accretion in the SN1987A data. 

Detecting Supernova Neutrinos Before describing the model, we recall how 
many events we expect from the most important detection reaction, the inverse beta 
decay (IBD): De +p -^ e^ +n. The number of events is: N^v — Np x Fj^^ x ap^p. The 
cross section for a De with average energy E = 15 MeV is CTp^p w G^pE^ « 10"'*"'^ cm^. 
The C-e fluence (time integrated flux) can be written by Fp^ = ''/^jji ^ — -^^ — r^ ~ 
10^° p^, where Eb — 3-10^'^ erg is the gravitational binding energy, D is the distance 
of the supernova and "6" are the neutrino types. For 1 kton detector, there are 
NpKil kton X 10^ X 6 ■ 10^^ x 2/18 x 10^^ protons. So the number of expected events 
is about 10. This rough estimation agrees with the number of events observed: 16 in 
Kamiokande-lP (2140 tons), 8 in IM^ (6800 tons) and 5 in Baksaiff^' (200 tons): 
29 events in 30 seconds, that include a few background events. 

Other reactions are expected to yield less events in water Cherenkov (as 
Kamiokande-II, 1MB, Super-Kamiokande) or scintillators (as Baksan, LVD, Kam- 
LAND). This is true, e.g., for the elastic scattering reaction, where the cross section 
(Tes ^ G^rrieE is much smaller since E ^ Trig. However the electrons are scattered 
by supernova neutrinos in the same direction of the incoming neutrinos, for the same 



^Other features of the new analysis: (i) energy, time and direction of each event are taken into ac- 
count; (ii) the correct background^J, finite detection efficiency and energy resolution are described; 
(iii) dead times and live time fraction are included; (iv) only the relative times are used and the 
delay of the detector response, called offset times, is accounted for; (v) frequentist techniques of 
inference are applied with an unbiased likelihoooii; and (vi) the full 30 s analysis window is consid- 
ered. An updated cross section of IBE^ is used, an improved description of the neutrino spectrum 
is introduced and neutrino oscillations are accounted for in a suitable approximation. 
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reason, E ^ mg. In other words, elastic scatterings will allow 'neutrino imaging' as 
is illustrated in the figure, where we show the directions of arrival of the supernova 
events. The cluster of events in the center, due to elastic scattering, is well visible 
over the background of inverse beta decay events, which are only mildly directional. 
The statistical analysis shows t hat the d irection of the supernova is determined with 
an accuracy of few degrees. Se c ^ ^ ^ for further discussion. 

We close this introductory note remarking that supernova neutrinos are a spe- 
cial target of neutrino astronomy. Indeed, despite the rarity of the observable core 
collapses, these neutrino can be certainly detected as demonstrated by SN1987A; 
moreover, they will permit to shed light on many open theoretical problems, re- 
garding several fields: astrophysics, nuclear physics and particle physics. 

Model of Neutrino Emission We describe now the model for the neutrino flux. 
Each emission phase is characterized by its intensity, its duration and the average 
energy of the emitted neutrinos. So, we have 6 astrophysical parameters, namely: 
the initial mass (Ma), the time scale (tq) and the initial temperature (Ta) for the 
accretion phase; the radius (Re), the time scale (tc) and the initial temperature 
(Tc) for the cooling phase. Moreover, in order to take into account the delay of 
the first observed event with respect to the first neutrino, we introduced other free 
parameters, called offset times. Since the clocks of Kamiokande-II, 1MB and Baksan 
were not synchronized, we need 3 offset times. For details, see^. 

The expectations on these astrophysical parameters are: Re ~ Rns = 10 — 20 km, 
Tc = 3 — 6 MeV and the duration of the cooling phase should be few (or many) 
seconds. The accretion phase has Ma < 0.6 Mq, Ta is the few MeV range and the 
accretion should last just ~ 0.5 s. The accretion Pg luminosity can be estimated by: 

Laccr - iV„(^e+n)r* ^ 5 X lO^^H!! ( ^I^] ( ^) ( ^^ , (1) 

n\ c+„/ a gg^ \0.1MqJ \Q.6j \2MeVj ' ^ ' 

where Nn is the number of neutrons, that can be expressed in term of an accretion 
mass. Ma, and of the fraction of neutrons in the environment, y„. The cooling i?e 
luminosity has instead the form: 

Leooi - R'X - 5 X lO^iH^ f-^]' (-^X , (2) 



sec 



lOkmJ V5MeVJ 



where Re and T^, are, respectively radius and temperature of the cooling. These 
formulae make evident that the two phases are described by very different physi- 
cal models, namely, by a transparent atmosphere and by an opaque (black body) 
radiator; this is the reason why the first phase is much more luminous. 

The two phases are not contemporaneous. We parameterize the D^ flux as follow: 

0Pe W =4'a + il~ Jk{t)) X (f>,{t - Ta). (3) 

Above jk represents a function that terminates the accretion phase at i ~ Ta, 
approximated by exp [— (i/ra)^]. The flux is dominated by the accretion phase at 
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t <^ Tai whereas the coohng phase begins at i ~ r^ and eventuaUy dominates the 
flux when t ^ Ta- With previous considerations in mind, we understand that Eq. ^ 
describes the passage from an emission of volume to an emission of surface. The 
spectrum is quasi-thermal at any time; furthermore, the neutrinos luminosity and 
their average energies are smooth functions of the time. 

The best fit values of the astrophysical parameters of accretion and cooling 
emission phases are found to be: 

R^ = W+l km Ta ^ 4.6to:6 MeV, r^ ^ 4.7ti:2 s, 

M, = 0.2210:68 Mo, Ta^2At"oiMeY, r, = O.SSt^';^? s. 

The large errors are due to the limited statistics. The results are close to what we 
expect from the standard collapse, in particular we get £b = 2.2 x 10^^ erg. There is 
a 2.5 a evidence for the accretion phase. The 11 early events of Kamiokande-II (6), 
1MB (3) and Baksan (2) have a great probability to belong to the accretion phaseO 

Gravitational Waves and Neutrinos Gravity Waves (GW) are predicted by 
general relativity. They have not been observed directly yet, but we will have soon 
detectors of enhanced sensitivity. Core collapse supernovae can emit GWs during 
the collapse (or during the explosion) of a core collapse SN due to the change 
of the quadrupole moment of the star structure. Recent simulations show that a 
gravitational signal is emitted when the collapse of the inner core halts, as dictated 
by the stiffening of the equation of state at nuclear density. The consequent bounce 
of the outer core is pressure dominated without strong influence of the rotation. 
Therefore, it is possible to define a generic GW waveform which exhibits a positive 
pre-bounce rise and a large negative peak, followed by a ring-down; so the time 
of the bounce is strongly correlated to the time of the maximum amplitude of the 
gravitational signal. 

The duration of the GW signal is about 10 ms. Therefore, to help the search of 
such signals, one would like to identify the time of the bounce with an er ror of the 
same order studying other types of signal emitted from this event. In ref.^LJJit was 
argued that it is possible to identify the time of the bounce within ^ 10 ms by an 
analysis of the Pe signal from the explosion of a galactic core collapse supernova; 
i.e., neutrinos can provide the required trigger for the search of GW. 

In fact, extensive simulations of core collapse SNe shows that the onset of i>e 
luminosity is closely related to the time of the bounce. The time of the bounce 
Tbounco can be determined by the following equation, where the times in uppercase 



^ We add a remark on backgrounds events. The 6"'-event of Kamiokande-II, with energy below 
7.5 MeV, has a probability of 85% to be due to background; a posteriori, it should not be attributed 
to accretion. Similarly, we found that the 13*'' — 16*'' events are almost surely due to background 
and there is still some chance of another background event. Similarly in Baksan, where the number 
of events is larger than expected and this is quite likely a priori. In absence of more precise 
information, we assumed that 1MB was background free; we checked that the inferences do not 
change significantly assuming that it had a background rate similar to the one of Kamiokande-II. 
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are absolute times, in UT, whereas those in lowercase are relative intervals of time: 

-t bounce ^ -i 1st (tew ~r imass ^ tfly I ^oftjj (.4j 

Tist is the time of the first neutrino event detected. The time icw is the interval 
between the bounce of the outer core on the inner core and the beginning of i?e 
emission. This is reliably known and ranges within <gw = 1-5 — 4.5 ms. The time 
imasa is thc delay, due to neutrino mass, between the arrival of GW and neutrino 
signal; if we impose the cosmological bound J^i ^^vi < 0.7 eV, this is negligible. The 
time interval ta^ is the time of fly between the two detectors and depends on the 
SN position in the sky. Finally the non-negative parameter t^^ is the difference in 
time between the first neutrino and the first event detected. In summary, the main 
terms in Eq. (j4]) are the fly time itiy and the offset (or response) time iotf- 

For the analysis of a future galactic supernova event, it is important to consider 
the finite rising time of the Pg signal. This can be done multiplying the flux of 
Eq. ([3]) by the function /r = 1 — e^*/'^''. The rising time r^. « 50 — 150 ms is a 
very important and new parameter of the astrophysical model. It is related with 
the initial production of v^ and depends strongly on the velocity of the shock wave. 

How to determine experimentally the time of fly from neutrinos? If we know 
astronomically the direction of the SN, it is easy to correct for the difference of 
arrival times. But even if we do not know it, we can rely on elastic scattering (ES) 
events, that are directional and suffice to determine the time of fly precisely enough. 
(Note in particular that this can be applied to a supernova without optical output). 

Thus, the problem reduces to the estimation of ioff, i-e., the delay of the response 
of the detector to the neutrino signal. If we have enough data and if we reconstruct 
at the same time t^, it is possible to reconstruct successfully t„ ,f by fitting the data 
to the expectations. As we can see from Table 3 of the reference^, these conditions 
are expected to be satisfied for a galactic supernova event: The response time and its 
error are correctly estimated by the analysis. We conclude that the future galactic 
supernova can provide us a precious information to test a key prediction of general 
relativity and note incidentally, again frorrM>JJ, that also the other astrophysical 
parameters will be reconstructed with very good precision. 

Neutrino Mass It is known since ZatsepirM^ that supernova neutrinos permit us 
to investigate the absolute mass scale of neutri nos. Sev eral works attempted this 
after SN1987A: see the Table for a partial revie w '•' II " i We will mostly discuss the 
last two entries of this Table, but before doing that, we recall the idea. 

The flux, Eq. ([3]) , is a function of the emission time ti , namely the time measured 
from the beginning of antincutrino emission. In terms of the absolute emission times 
we have ti — T? — Tq , that can be rewritten taking into account the velocity of the 
neutrino Vi and the absolute detection time T^ as follows: 
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Reference m^ CL 



Nat 326, 476 


< 11 eV 


- 


Nat 329, 689 


0(13) eV 


- 


PRL 58, 1906 


< 12 eV 


- 


PRL 58, 2722 


< 26 eV 


- 


PRL 59, 1864 


<5eV 


- 


PLB 200, 366 


< 16 eV 


95% 


PRD 35, 3598 


< 13.5 eV 


la 


ELett 4, 953 


< 5.7 eV 


- 


MPLA 2, 905 


3.4 ±0.6 eV 


- 


CNPP 17, 239 


< 30 eV 


- 


BIHEP-CR-87-01 


4.5±0.9eV 


- 


DTP/87/12 


3.4 ±0.5 eV 


- 


NPB 299, 734 


m^ = 4tf3 eV2 


- 


PLB 196, 259 


< 10 eV 


- 


PRD 41, 682 


< 14 eV 


95% 


NPB 437, 243 


< 19.6 eV 


95% 


PRD 65, 063002 


< 5.7 eV 


95% 


iarXiv:1002.3349| 


< 5.8 eV 


95% 



where T^ is the minimum, possible detection time. The first term in the r.h.s. 
T/ — Ti = 5ti is known from the experiment. The second one, T^ — Tq = i„ff is by 
definition the offset time, namely the delay between the first observed event and 
the moment when the first neutrino has possibly reached the detector. The last 
one, denoted by Ati, describes the effect of the finite neutrino mass. Its numerical 
expression is: 

D /jTz^c^ \ VlO McV ^^ 



At, = 2.6 sec X -^ — . (6) 

50kpc VlOeVy V E^.^ J 

Putting together the above definitions, the emission time of each event, that enters 
the likelihood through the antineutrino flux, will be written as: 

ti = SU + t^::~AU. (7) 

Again, the first term in the r.h.s. is the relative time between the i-th and the first 
observed event in the considered detector, known directly from the data without 
significant error. The last two terms, instead, have to be estimated by fitting the 
data; both of them are positive (and thus lead to some cancellation) but depend in 
a different way from the neutrino energy, see Eq. (j6]). E^^^i, in turn, can be inferred 
from the measured energy of t he p ositron, Ei, which is known up to its error 6Ei. 
From our statistical analysis^^, we obtain from SN1987A data the bound 

m^ < 5.8 eV at 95% CL. (8) 

The existence of a phase of accretion of about 0.5 s explains why we are able 
to probe such neutrino masses. From Eq. ([6]) it should be evident that the most 
important information to determine the mass is contained in the first, low energy 
events. This consideration selects as most relevant the events of Kamiokande-II, 
which incidentally, are also the events that probe the existence of the accretion 
phase. The astrophysical uncertainties are not very relevant; instead, the fact that 
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offset time is unknown contribut es to worsen the bound. Note that the agreement 
of the last two entries of the tabi c I 4s to a large extent accidental; e.g., the bound 
changes by 5 — 10% adopting other conventional statistical procedures. 

By mean of simulated data, we demonstrated the possibility to probe sub-eV 
neutrino masses with Super-Kamiokande after a future galactic supernova. As dis- 
cussed iii^, this limit requires a very precise knowledge of the time and it is subject 
to strong fluctuations, related to the position of the first low energy event. 

Remark on the Relevant Time Scales We would like to summarize here the 
above results emphasizing the role of the relevant 'time scales': 

1) We recalled that the conventional astrophysical picture of the explosion contains 
a relatively short time scale, namely the duration of the accretion Ta ^ 0.5 s. We ar- 
gued that SN1987A data already provide some evidence for this time scale, thanks 
to the fact that almost 40% of the observed events fall in this phaseO A future 
galactic supernova will give us a much better determination of the astrophysics and 
plausibly also of 10 ms (or even few ms) structures in the signal. 

2) The search for GW can profit of the detection of neutrinos. We remarked that 
the observation of a future galactic supernova neutrino event could permit a deter- 
mination with 10 ms precision of the moment of the bounce, when a burst of GW is 
plausibly emitted. (The same cannot be done for SN1987A, since the absolute time 
of the events is unknown, except for 1MB; thus, the bounce should have occurred 
0.76 s before the first event seen by 1MB, see Eq. (32) of ). The postulated rise 
time of the signal, of the order of 50-150 ms, can also be measured. 

3) The determination of the neutrino mass has to go through a precise determina- 
tion of times: see Eq. ([7]); the smaller the time scale we probe, the better the limit 
on the mass we obtain. E.g., the bound from SN1987A (several eV) is essentially 
determined by the existence of an accretion phase, while for a future supernova the 
relevant time scale will be the rise time of the signal (probing the sub-eV region). 
This could be further improved detecting a hypothetical shorter burst of Pg, lasting 
only few ms. This is not expected to exist for a standard collapse, but it has been 
found in the first numerical simulations concerning the formation of a hybrid stap^. 

Conclusions The problem to explain SN explosion is still open, however, a ref- 
erence (standard) model does exist. In this work, we focussed on this standard 
scenario and discussed possible observational tests. 

The SN1987A data present unexpected features, however, KII, 1MB and Baksan 
data fit in a specific model. They show a hint of an initial high-luminosity phase, as 
the one expected for a standard neutrino emission. A future galactic SN will permit 
much more precise tests providing a huge amount of new information; remarkably. 



'^This is a significant extension of the usual approach to SN1987A data analysis (as summarized, 
e.g., in Bahcall book) where the neutrinos are thought to originate from a smooth, thermal emission 
with a much longer time scale, Tc = several seconds. 
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it will also give the GW burst timing within ^ 10 msec. Also it will be possible to 
establish a relatively tight neutrino mass bound, up to the sub-eV range. 

Finally, we noted that many of the interesting results on astrophysics and particle 
physics that can be obtained by observations of supernova neutrinos are essentially 
based on precise measurements of time. This consideration emphasizes the crucial 
importance to improve our knowledge on the astrophysics of the neutrino emission. 
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